ABSTRACT 



I study phase transitions occuring in noncollinear magnets by means of a 
self-consistent screening approximation. The Ginzburg-Landau theory involves 
two N-component vector fields with two independent quartic couplings allowing 
a symmetry-breaking scheme which is SO(N) x SO (2) -> SO(N-2) x SO(2) diag . I 
find that there is a second-order phase transition in the physical cases N=2,3, D=3 
and that there is no fluctuation-induced first-order transition. This is very similar 
to the case of the normal-to-superconducting phase transition as recently found by 
Radzihovsky. The exponents are r](N = 3, D = 3) as 0.11, r](N = 2, D = 3) as 0.15 
and go smoothly to the large-N limit. 
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Many magnetic systems have a low-temperature ordered phase that breaks completely 
the rotation invariance due to a noncollinear pattern of the magnetic moments 1 . Exam- 
ples includes the rare-earths Ho, Dy and Tb. This is generic for vector-spin models with 
strong enough competing interactions. The phase transition associated with this order- 
ing does not belong a priori to the well-studied universality classes corresponding to the 
collinear ordering with symmetry breaking O(N) — > O(N-l). Several groups have investi- 
gated the critical behaviour of the three-dimensional stacked triangular (STA) Heisenberg 
antiferromagnet which is a simple example of commensurate noncollinear ordering and 
there is general agreement 2 ' 3 ' 4 that this system has a second-order phase transition with 
exponents 3 v = 0.585(9), r y/u= 2.011(14) that do not correspond to the O(N) exponents. 
This peculiar set of exponents is apparently associated with a new universality class: they 
appear also in the body-centered tetragonal antiferromagnet 5 . This class should also in- 
clude the transition 1 from the 3 He liquid to the A-phase, Josephson-j unction arrays in a 
transverse field as well as the fully frustrated bipartite lattice (Villain model) . In the XY 
case there is also a new set of exponents 1 ' 2 : v = 0.54(2), 7 = 1.13(5) 

The renormalization group has been applied to the Ginzburg-Landau theory of such 
systems by Garel and Pfeuty 6 . They used the expansion in e = 4 — D, where D is the 
dimension, for any number N of components of the vector spins. In this framework, there is 
a Heisenberg 0(2N) fixed point which is always unstable. In the physical case N=3 there 
is no new fixed point and this runaway is interpreted as indicating a fluctuation-induced 
first-order phase transition. This is curiously similar to the normal-superconducting (NS) 
phase transition as revealed by Halperin, Lubensky and Ma 7 . 

Recently, the NS phase transition has been studied by Radzihovsky 8 who used the 
self-consistent screening approximation (SCSA) due to Bray 9 . This approximation includes 
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an infinite subset of the 1 /N expansion by means of a simple self-consistency condition on 
the propagator. It allows to obtain the exponent rj as a function of N and D. In the 
NS case the SCSA leads to a stable fixed point in the physical case D=3, in agreement 
with Monte-Carlo findings and other theoretical arguments 10 . In the neighborhood of 
D=4, the corresponding fixed point survives for all values of the number of components 
of the order parameter, contrary to the prediction of the 4 — e expansion, suggesting that 
the fluctuation-induced first-order transition is an artifact of the e expansion. It is the 
purpose of this Letter to apply this very same method to the Ginzburg-Landau theory of 
helimagnets. I find that there is a non-trivial fixed point that survives in the whole domain 
2 < D < 4 and 2 < N < oo which is different from the 0(2N) Heisenberg fixed point. In 
the physical cases of interest, I find r](N = 3, D = 3) « 0.11 and i](N = 2, D = 3) 0.15. 
This fixed point goes smoothly to the N = oo result and there is no fluctuation-induced 
first-order phase transition. This is very close to the analysis of the NS transition of Ref.8. 
The mere existence of this fixed point for N = 2, 3, D = 3 provides a natural explanation 
to the Monte-Carlo results and, possibly, of some experimental results. 

The Ginzburg-Landau theory for a generic Heisenberg helimagnet involves two vec- 
tor fields that correspond to the Fourier modes of the magnetization near the ordering 
wavevectors ±Q. The effective action contains two quartic invariants: 

A= \ (( V ^) 2 + ( V ^ 2 ) 2 ) +r + + 

+U (jfi + + V ((<?! • ^ 2 ) 2 " ■ (1) 

This free energy has a global symmetry O(N) x 0(2). When the coefficient v is positive, the 
ground state consists of orthogonal vectors and the residual symmetry is 0(N— 2) xO(2)^ a£ , 
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where the subgroup 0(2) diag acts diagonally on vector indices and internal (1,2) indices 11 . 
The global symmetry allows only the two quartic invariants in Eq.(l). A detailed study of 
this model has been performed by H. Kawamura 12 . In the D = 4 — e calculation, one finds 
the 0(2N) Wilson-Fisher fixed point on the line v = at a distance e from the origin, in 
the (u,v) plane. The operator ~~ ($i ' <p2^j ] opens a direction of instability for all 

values of N. However, if N is greater than N C (D) = 21.8 — 23. 4e + 0(e 2 ) (obtained from 
a two-loop computation 12 ), there is an additional stable fixed point with u* ^ 0, v* ^ 0. 
In the neighborhood of the upper critical dimension, there is thus a dividing line N C (D) 
in the (N, D) plane above which one has a second-order phase transition and below which 
one conjectures a fluctuation-induced first-order phase transition. This is similar to the 
normal to superconductor transition 7 ' 8 except for numerical factors: there the transition 
is second-order only when the number of complex components of the order parameter is 
larger than 183. 

The self-consistent screening approximation is an improvement over the large-N expan- 
sion which has been applied with success to several physical problems 8 ' 9 ' 13 . In the leading 
large-N expression for the self-energy of the basic field containing the geometric sum of 
bubbles, one uses renormalized propagators everywhere: this leads to a self-consistency 
equation that goes beyond the simple result N — > oo. An interesting byproduct is that 
this approximation gives back the N — > oo result automatically. Its validity extends to all 
N and D values but, of course, this is not a systematic procedure and we do not expect 
to get precise numerical values for the exponent rj. The propagator of the theory (1) is 



defined as Gn(fc) = ((f>i{k)(f>f{-k)) = (0£(fc)0£(-AO) = G 22 (k) = G(k). The self-energy 
is defined by G~ 1 (k) = k 2 + r + E(fc). The Dyson equation is then: 
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where the free propagator is G 1 (k) = k 2 + r and the dressed vertices ii(p), v(p) are given 
by series involving only powers of the polarization bubble H(k) = J G(p)G(p + k): 

m - m= u + v(2u-v/2)NU(p) 

yP) l + vNTl(p)' yP) l + 4uNU(p) + v(4u-v)N 2 U(p) 2 ' {) 

Note that it is the full propagators that enter the quantities n(p), u(p), v(p) and thus Eq.(2) 
is a self-consistency condition. The corresponding equation in the large-N limit is obtained 
by using the free propagator Go in the right-hand side of Eq.(2). 

The SCSA strategy 9 amounts to considering criticality i.e. r = and writing the 
scaling form of the propagator G~ 1 (k) = k^k 2-11 since there is no longer any correlation 
length. With this scaling form for G one can compute the polarization bubble following 
Bray 9 : 

ll[P) ~ (4n)»/ 2 T(D-2 + V )T 2 (l- V /2) * ^ V ' W 

The Heisenberg 0(2N) fixed point can be found by setting v = 0. If we suppose that 
IT blows up for small momenta then u(p) m l/4A^n. One can now match the powers of 
k~ 2v in Eq.(2): this fixes the value of rj through the condition: 

N= IW2 - l)r(2 - y)T(D - 2 + 77)1X1 - 77/2) 

" T(D/2 + 7]- 2)T(D/2 + V /2 - l)T(D/2 - V /2 + 1)T(2 - 77 - D/2) ' 1 ' 

This is precisely the result of the SCSA for the 0(2iV)-vector model as expected. As long 
as D < 4, one can check the assumption that 77 leads to IT >> 1 at low momenta (when 
D > 4, then IT is negligible and one encounters only the Gaussian fixed point). 

The other fixed point corresponds to u 7^ 0: assuming also IT >> 1 we use the 
asymptotic behaviour v(p) ~ 1/iVTI and u(p) ~ 1/2./VTI, deduced from the series (3). 
These leading terms add up to give a factor of 3/iV in the right-hand side of the Dyson 
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equation and thus the only modification with respect to the pure Heisenberg case is that 
in Eq(5) the factor N should be replaced by N/3: 

N = T(y/2 - l)r(2 - y)T(D - 2 + V )T(1 - V /2) 

3 T(D/2 + i 1 -2)T(D/2 + i 1 /2-l)T(D/2-r 1 /2 + l)T(2-r 1 - D/2)' 1 ' 

This equation leads to a rjscsA which is well-behaved in the physical case N=2,3, D=3. 
It can be solved now numerically or expanded in various limits. Numerically t]scsa(N = 
3,D = 3) ~ 0.11 and t)scsa(N = 2,D = 3) ~ 0.15, a perfectly sensible result. Present 
Monte-Carlo estimates 3 favor a smaller r] for N=3 since 7/z/ = 2.011(14). The SCSA 
should not be expected to be quantitative: this is the mere existence of the fixed point 
which is the relevant information. 

The large-N limit of the SCSA agrees by construction with the result of the direct large- 
N result 12 : one recovers the result i] = 6[4/D — 1]Sd/N where Sd = sm[jr(D — 2)/2]T(D — 
l)/(2nT 2 (D/2)). In the neighborhood of D=4, one finds vscsa = 3e 2 /4iV + 0(e 3 ). This 
e-expansion is well-behaved for any N: in the SCSA there is no hint of a fluctuation-induced 
first-order transition. This limiting case matches the e-expansion only if one takes also the 
limit iV — > 00. In fact, r]4- e = f(N)e 2 + 0(e 3 ) with f(N) a complicated function of N 
which has the asymptotic behaviour f(N) ss 3/4A^ as — > 00. 

In the case of the NS transition, Radzihovsky has shown 8 that rises A is well-behaved 
near D=4 but has a singularity in the e-expansion. This has led him to propose that 
the lack of stable fixed point seen in the e-expansion studies 7 is due to a breakdown of 
the expansion itself instead of a fluctuation-induced first-order phase transition. Since it 
is difficult to estimate the validity of the SCSA itself, it may be that the SCSA is not 
accurate enough to capture the fluctuation-induced transition. In the present context of 
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noncollinear magnets, the situation is slightly different: rises A is well-behaved near D=4 
and has a regular e-expansion for any N. 

In the neighborhood of the lower critical dimension, the exponent rj can be expanded 
in D=2+e with the result 7]scsa = 3e/2(A — 3) + 0(e 2 ). There is a singularity at N=3 
in the e-expansion due to the fact that r\ no longer vanishes (in the SCSA) at D=2 when 
A < 3. For N=3, one has from the SCSA (equation 6) 1]=0 at D=2 which is coherent with 
the fact that there is no phase transition. This behaviour appears of course in the SCSA 
treatment of the 0(N)-vector model since Eqs 5 and 6 differ by the substitution A — > A/3. 
However in the O(N) case, rises A = e/(A — 2) near D=2 and rises A is nonzero at D=2 
when A < 2. Since we expect in a N-vector model that in D=2 r\ vanishes for A > 2, this 
means that the N-dependence from SCSA cannot be blindly believed near N=3. 

Noncollinear magnets have been studied near two dimensions in a sigma model 
approach 11 . The symmetry breaking pattern defines a homogeneous non-symmetric man- 
ifold G/H = O(A) x 0(2)/0(N — 2) x 0(2)di ag that specifies uniquely a nonlinear sigma 
model. For any finite N there is a fixed point in D=2+e which merges smoothly with the 
fixed point obtained in the large-N limit of the linear theory. The exponent r\ is given 
by 14 ' 11 : 

3A 2 -10A + 9 2 . 
^= 2(A-2)3 e + ^ 

For large N, this agrees with rises A ~ 3e/2A. In addition the formula (7) has the correct 

feature that it blows up for N=2 instead of N=3 as found in the SCSA. For N=3, the sigma 

model has a peculiarity: since 0(3) x 0(3) = 0(4), it describes 0(4) critical behaviour. 

This is incompatible with numerical simulations 3 in D=3. The reason for this failure is 

not known at the present time. A view has been reported 14 that the inability of the 
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e = D — 2 expansion to detect the nontrivial topological structure of the order-parameter 
space might be the reason for this failure. It is clear that the SCSA cannot shed any light 
on this problem since its functional dependence upon N is too approximate. In the SCSA 
there is no hint of any symmetry enhancement since 0(4) critical behaviour cannot be 
reached by the linear theory (1). 

In this Letter, I have studied the phase transition occuring in noncollinear magnets 
that break the full rotation group at low temperatures. A self-consistent screening approx- 
imation valid for all dimensions and any number of components of the order parameter 
leads to a second-order phase transition in the physical case of Heisenberg and XY systems 
in D=3. There is no sign of fluctuation-induced first-order phase transition contrary to the 
result of the standard 4 — e-expansion. This is in agreement with Monte-Carlo results on 
model systems that possess the correct symmetry-breaking pattern. There are also some 
experimental results 11 that may be explained by the existence of the corresponding uni- 
versality class. It is worth pointing out that the SCSA approximation is not a controlled 
approximation with a systematic expansion parameter and thus it is difficult to estimate 
its validity. The whole picture is close to the normal-to-superconducting phase transition 
as studied by Radzihovsky: here the SCSA leads also to a second order transition in the 
whole (N, D)-plane instead of the first-order transition predicted by the 4 — e-expansion. 
There is however a difference: in the NS case i]scsa is well-behaved but has a singularity 
when expanded in e while in helical magnets even the expansion is regular. In the two 
cases (NS and helical magnets) the transition is continuous near two dimensions 11 ' 15 and 
in D=3: this is reproduced by the SCSA calculation. 

Acknowledgments: I thank F. David, E. Guitter and J. Zinn- Justin for discussions. 
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